Hyperfine structure of the 2 3 5'— state in the 9 Be atom. 
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Abstract 

A number of bound states properties and the hyperfine structure of the 2 3 S(L = 0)— state in 
the Be atom are determined from the results of highly accurate bound state computations. The 
energies of the hyperfine structure levels for this state are e(F = I ) = -13725.612 MHz, e(F = |) 
= -5490.245 MHz and e(F = f ) = 8235.367 MHz, respectively. The hyperfine structure splittings 
are A 12 = 8235.367 MHz and A 23 = 13725.612 MHz, respectively. 



*E-mail address: afrolov@uwo.ca 



1 



In this study we investigate the hyperfine structure of the triplet bound S(L = 0)— states 
in the four-electron Be atom(s). As follows from the results of astrophysical observations 
the triplet states in the four-electron 9 Be atom manifest an interesting hyperfine structure. 
The main goal of this study is to investigate this hyperfine structure and determine the 
corresponding hyperfine structure splittings which can be measured to very high accuracy 
in modern atomic experiments. Here we restrict our analysis to the lowest triplet 3 S(L = 
0)— state in some Be atoms. It should be mentioned that the bound triplet states in different 
four-electron atoms and ions have been neglected for quite some time. Various theoretical, 
computational and experimental works performed recently for the four-electron atoms and 
ions were mainly oriented to the singlet states in such systems. For the triplet bound states 
only a few highly accurate results are now known (see, e.g., and references therein). 

On the other hand, different triplet (bound) states in the four-electron atoms and ions are 
of great interest in a large number of physical problems, including spectral analysis of the 
stellar and laboratory plasmas, accurate prediction of the properties of the light element 
plasmas at high temperatures and arbitrary pressures, etc. 

The Hamiltonian H of the four-electron atomic system is written in the form 
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where h is the reduced Planck constant, m e is the electron mass, e is the absolute value of 
the electron charge of the electron. Also, in this equation Q and M n are the electric charge 
and mass of the nucleus expressed in e and m e , respectively. It is clear that M n 3> 1 and 
to form some bound states in the four-electron atoms/ions the following condition must be 
obeyed Q > 3. Here and everywhere below the subscript 5 denotes the atomic nucleus, while 
subscripts 1, 2, 3 and 4 stand for electrons. The four-electron Be atom has two independent 
series of bound states: singlet states and triplet states. The ground singlet state with L = 
(also called the l^S^L = 0)— state) has the minimal total energy. The total energies of the 
two bound P(L = 1)— states (singlet and triplet) are well above the energy of the ground 
l 1 ^— state (see, e.g., [3j). In turn, the bound 2 3 S— state of the Be atom is located above 
these two P(L = l)-states. The electron configuration of this state is ls 2 2s3s. Note also 
that the total energy of bound states in the four-electron Be atom must be lower than the 
total energy of the ground (doublet) 2 S— state of the three-electron Be + ion which equals 
~-14. 324763 17647 a.u. Otherwise, any discussion of these states is not correct. Our best 



total energy determined for the bound 2 3 S— state of the °°Be atom is -14.43005958 a.u. To 
obtain this energy we have used up to 1600 basis functions (see below) in our calculations. 

The wave function of the bound 2 3 S(L = 0)— state in the Be atom is represented in the 
following form 

* = A e [ip(A; {rij})(a/3 - 0a)aa] + A e [<f){B; {r i:j })(2aa(3 - (3aa - a(3a)a>] (2) 

where the notation {r^} designates all ten relative coordinates (electron-nuclear and 
electron-electron coordinates), while the notation A e means the complete four-electron sym- 



metrizer, while notations a and (3 denote spin- up and spin-down functions ^|. The explicit 
formula for the A e operator is 

A e = e — P\2 — P 13 — p23 — P\A — P24 — P34 + P 123 + f 132 + P.24 + Pl42 + P 134 + Pl43 

+-P234 + P243 — -P1234 — -P1243 ~ -P1324 — -P1342 — -P1423 ~ -P1432 + -P12-P34 + P13P24 + -P14-P23 (3) 

Here e is the identity permutation, while Py is the permutation of the zth and jth identical 
particles. Analogously, the notations and Pijki stand for permutations of the three 
and four identical particles (electrons). In real calculations one needs to know the explicit 
expressions for the spatial projectors only. These projectors are usualy obtained by applying 
the A e operator to each component of the wave function in Eq.(J2]). At the second step 
we need to determine the scalar product of the result and incident spin function. After 
the integration over all spin variables one finds the formula for the spatial projector. For 
instance, in the case of the first term in Eq.(J2]) we obtain 

'Pipip — ^ /g(2e + 2P 12 — P 13 — P23 — Pu — P24 — 2P 34 — 2P 12 P 34 — P 12 3 — A24 — A32 

— P142 + P134 + -P143 + -P234 + -P243 + -P1234 + -P1243 + -P342 + P1432) (4) 



The functions ip(A; {rij}) and <p{B] {rjj}) in Eq.© are the radial parts (or components) of 
the total wave function For the bound states in various five-body systems these functions 
are approximated with the use of the KT-variational expansion written in ten-dimensional 



gaussoids j|3], e.g., for the ip(A; {r^}) function 
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where Na is the total number of terms, Ck are the linear variational coefficients and V is 
the spatial projector defined in Eq.(j4"]). The notations A and B in Eq.([2]) stand for the cor- 
responding set of the non-linear parameters {a^} and {bff} in the radial wave functions, 



Eq.(j5]). It is assumed here that these two sets of non-linear parameters are optimized inde- 
pendently of each other. In general, the KT- variational expansion is very effective for various 
few-body systems known in atomic, molecular and nuclear physics. A large number of fast 
algorithms developed recently for optimization of the non-linear parameters in the trial wave 
functions, Eq.flS]), allow one to approximate the total energies and variational wave functions 
to high accuracy. In some cases, however, the overall convergence rate for some bound state 
properties is sufficiently lower than for the total energies and other properties. In particular, 
it was found that the expectation values of the two- and three-particle delta-functions, e.g., 
(5(r e jv)), converge very slowly. In other words, it takes a long time to approximate these 
expectation values to high accuracy. The total energy and other expectation values do not 
change drastically during such an additional optimization of the non-linear parameters in 
the wave functions. In this study we report a number of improved expectation values for the 
bound triplet 2 3 S— state of the Be atom (see Table I). The improved expectation value of 
the electron- nuclear delta-function (5(r e jv)) from Table I allows one to make more accurate 
predictions of the hyperfme structure splitting for this state in the 9 Be atom (see below). 

The total electron spin of the triplet 2 3 S— state of the Be atom is different from zero. 
This means a possibility to observe a spin-spin interaction between the total electron spin 
S e and nuclear spin I/y in those atoms where | I n |^ 0. For instance, for the triplet state in 
the 9 Be atom we have S e = 1, Ijy = | and we can see the triplet of hyperfme structure states 
with F = §,§,§, where F is the absolute value of the total angular momentum operator 
F = L + S = S e + lAr for the whole atom, i.e. the sum of the electron and nuclear spins 
(for the triplet S— tates). In the first approximation, the levels of hyperfme structure are 
represented by the Fermi- Segre formula 

Hhf = ^^ B ^N9e9N(S(r eN ))^[F(F + 1) - I N (I N + 1) - S e (S e + 1)] (6) 

where a = 7.297352568 • is the dimensionless fine structure constant, hb = \ is the 
Bohr magneton (in atomic units) and fi^ — ^b^-, where ^ = 1836.15267261 is the ratio of 
;he proton and electron masses. The electron gyromagnetic ratio g e = -2.0022 319304 3718 
6]. The factor gjy for the 9 Be nucleus is -1.17761 (|) = -0.78507. Finally, the formula for 
the Hhf takes the form 

Hhf (MHz) = 314.121190927 ■ (8(r eN )) ■ [F(F + 1) - I N (I N + 1) - S e (S e + 1)] (7) 
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where we used the factor 6.57968392061-10 9 (MHz/a.u.) to recalculate the result from the 
atomic units to MegaHertz. By using the expectation value of the electron-nuclear delta 
function (5(r eN )) sa 8.739055 a.u. one finds the following formula for H h f(MHz) 

H h f(MHz) = 2745.12237[F(F + 1) - I N (I N + 1) - S e (S e + 1)] = 2745.12237 x 

[F(F + 1) - f ] (8) 

since for the triplet 2 3 S'— state in the 9 Be atom we have S e = 1 and Jjv = §■ 

The energies of the hyperfme structure levels are e(F — |) = -13725.612 MHz, e(F = |) 
= -5490.245 MHz and e(F = |) = 8235.367 MHz, respectively. The corresponding hyper- 
fine structure splittings are A 12 — 8235.367 MHz and A23 = 13725.612 MHz, respectively. 
These values are significanly more accurate than our hyperfme splittings reported in [l]. 
Note also that these hyperfme hyperfme structure splittings can be measured in the modern 
experiments. Analogous problem for the 2 3 P(L = 1)— state of the 9 Be atom has never been 
considered with the use of highly accurate variational wave functions. 
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TABLE I: The expectation values of some electron-nuclear (en) and electron-electron (ee) proper- 
ties of the 2 3 S— state of the °°Be atom (in atomic units). 



(rj) 




(reN) (r 2 eN ) (r% N ) (rj N ) 


14.274856 


2.035997 


2.632627 17.17728 148.3102 1468.243 




(ree 1 ) 


(ree) (r 2 ee ) (rl) (r e 4 e ) 


1.503444 


0.619276 


4.712471 35.19267 319.9908 3276.500 


(bl) 


(\p 2 n) 


($eN) (See) ($eeN) (^eeeN) 


3.6075710 14.8929525 


8.739055 0.39752 38.1614 0.0 
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